The dynamics of thin vibrated granular layers 
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We describe a series of experiments and computer simulations on vibrated granular media in a 
geometry chosen to eliminate gravitationally induced settling. The system consists of a collection 
of identical spherical particles on a horizontal plate vibrating vertically, with or without a confining 
lid. Previously reported results are reviewed, including the observation of homogeneous, disordered 
liquid-like states, an instability to a 'collapse' of motionless spheres on a perfect hexagonal lattice, 
and a fluctuating, hexagonally ordered state. In the presence of a confining lid we see a variety 
of solid phases at high densities and relatively high vibration amplitudes, several of which are 
reported for the flrst time in this article. The phase behavior of the system is closely related to that 
observed in confined hard-sphere colloidal suspensions in equilibrium, but with modifications due 
to the effects of the forcing and dissipation. We also review measurements of velocity distributions, 
which range from Maxwellian to strongly non-Maxwellian depending on the experimental parameter 
values. We describe measurements of spatial velocity correlations that show a clear dependence on 
the mechanism of energy injection. We also report new measurements of the velocity autocorrelation 
function in the granular layer and show that increased inelasticity leads to enhanced particle self- 
diffusion. 

PACS numbers: 45.70.-n,05.70.Fh,05.70.Ln 



INTRODUCTION 



Despite its apparent simplicity, a granular medium 
consisting of a collection of dry, cohesionless, identical 
spherical particles can display a wide range of complex 
behaviors for which no reliable mathematical models ex- 
ist. The empirical relations used by engineers to de- 
scribe granular flows cannot be reliably extended to deal 
with new processes or new environments. Attempts to 
improve models of granular media have been thwarted 
by a range of challenging problems. Collisions between 
granular particles may dissipate a considerable fraction 
of the available kinetic energy, so the results of equi- 
librium kinetic theory are not directly applicable. Fur- 
thermore, dissipation often leads to the appearance of 
large gradients and shocks, which makes it difficult to 
employ continuum equations 0, Another set of com- 
plications arises from the fact that macroscopic particles 
settle rapidly under the Earth's gravity, so that granular 
gases, in which particles interact through brief collisions 
rather than enduring contacts, are difficult to produce ex- 
perimentally. When, as is most often the case, collisional 
pressures aren't sufficient to support the weight of the 
granular media, a complex network of particle contacts 
results. This leads to additional complications because 
of the complexity of the frictional inter-particle contacts. 
Moreover, the inability of particles to surmount the grav- 
itational barriers preventing one particle from moving 
over another results in a kind of frozen-in metastable 
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state. This state has much in common with other glassy 
states of matter, including the fact that it is not very well 
understood. 

In this article we describe a series of experiments and 
computer simulations on vibrated granular media in a ge- 
ometry chosen to eliminate gravitationally induced set- 
tling. The system consists of a collection of identical 
spherical particles on a horizontal plate vibrating verti- 
cally, with or without a confining lid (Fig^. This sys- 
tem is well suited to investigate many of the complex 
non-equilibrium effects observed in excited granular me- 
dia and is simple enough that it can be accurately mod- 
elled in molecular dynamics simulations. By controlling 
the amplitude and frequency of the plate vibration, the 
density of spheres on the plate, the surface properties of 
the bottom plate, the size of the gap between the bot- 
tom plate and the confining lid, and the inelasticity of 
the spheres (by using different metals), we can observe 
a wide range of behaviors. At low densities and mod- 
erate vibration amplitudes, a horizontally homogeneous, 
disordered liquid-like state is observed. As the shaking 
amplitude is decreased, the granular liquid becomes un- 
stable to a 'collapse' of motionless spheres into a perfect 
hexagonal lattice. As described in Section „ this transi- 
tion is probably due to a bi-stability in the interaction 
between the balls and the vibrating plate at low vibration 
amplitudes. At higher densities we see numerous differ- 
ent solid phases, several of which are reported for the 
first time in this article. The phase behavior of the sys- 
tem is closely related to that observed in confined hard- 
sphere colloidal suspensions in equilibrium. The equilib- 
rium transitions are driven by entropy maximization, and 
the close similarity strongly suggests a common mecha- 
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nism. However we also observe nonequilibrium effects, 
including non-equipartition of kinetic energy between co- 
existing phases H, that directly affect the phase diagram. 

We have observed that the liquid-like states can have 
strongly non-Maxwellian velocity distributions, while in 
other conditions the distributions are indistinguishable 
from a Maxwell distribution0,0,0,0|- The granular gas 
displays strong spatial velocity correlations, but those 
correlations change dramatically when the smooth hori- 
zontal plate is replaced with a rough one0- Taken to- 
gether these results demonstrate that the statistical prop- 
erties of the granular gas are dependent on the external 
forcing that is necessary to maintain a steady state in 
the presence of dissipative interparticle collisions. This 
effect is probably exaggerated in the very thin granular 
layers that we have studied because every particle fre- 
quently interacts with the confining boundaries. Other 
experiments and computer simulations on quasi two- 
dimensional granular gasses have shown strongly non- 
Maxwellian velocity distributions and correlations that 
can be understood at least in part through consideration 
of the interactions with the boundaries |4 Isllsll^ Hol Ull 
This is consistent with theoretical work on 
the non-equilibrium steady state obtained when energy is 
supplied to the granular ga s by a spatially homogeneous 
random external forcing |l6lll7lll8lll9.] . Non-Maxwellian 
velocity distributions and algebraically decaying velocity 
correlations arise as a direct consequence of the energy 
injection. These distributions and correlations are neces- 
sary ingredients of a kinetic theory of granular fluids, and 
can have significant effects on transport coefficients. As 
an example, in this paper we report new measurements 
of the velocity autocorrelation function in the granular 
layer, and show that increased inelasticity leads to in- 
creased correlations, which in turn enhances particle self- 
diffusion. 

This paper is organized as follows: Sectioun (Methods) 
describes the experimental apparatus and measurement 
techniques as well as the model used for the molecular 
dynamics simulations. Section □ describes the phases and 
phase transitions observed. Section g reviews the statis- 
tical characterization of granular gasses. Finally, a few 
concluding remarks are offered in Section „ 



0.5%. For all of the experiments described in this paper 
the particles are smooth spheres of either stainless steel 
(coefficient of restitution e = 0.95) or brass (e = 0.77,) 
with diameters, cr, ranging from 1.5 mm to 3 mm, de- 
pending on the experiment. Most experiments and sim- 
ulations have been performed with a smooth plate, but 
in some cases we have used a 'rough plate', formed by 
a hexagonally close-packed layer of 1.19 mm diameter 
stainless steel spheres glued to the bottom plate. The 
particles are surrounded by an aluminum rim, and for 
some experiments they are confined from above with a 
thick, anti-static coated Plexiglas lid. The particles are 
illuminated from above to produce a small bright spot 
at the top of each of them when viewed through a video 
camera mounted directly above the plate. For all of the 
experiments described below, the plate vibration is sinu- 
soidal, and typically in the range oi ly — 30 — 100 Hz. 
The peak plate acceleration, Aui'^, where A is the am- 
plitude and Lu = 2'Ki', typically ranges from 1 gi up to 5 
g. It is often convenient to report results in terms of the 
dimensionless parameter F = Auj'^/g, 
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FIG. 1: Sketch of the experimental apparatus used in these 
studies. In some experiments, a rough bottom plate is used 
in place of the smooth one pictured here and in others, the 
confining lid is absent. In section □ the effect of varying the 
gap height, H, is described. 



METHODS 
Experimental Apparatus 

The experimental apparatus consists of a rigid alu- 
minum plate approximately 20 cm in diameter that is 
mounted horizontally on an electromagnetic shaker that 
oscillates the plate vertically. The plate is carefully lev- 
elled, and the peak acceleration, monitored with a fast- 
response accelerometer mounted on the bottom surface 
of the plate is uniform across the plate to better that 



Several digital video cameras have been used for data 
acquisition, including a high-resolution camera for study- 
ing spatial correlations (Pulnix-TM1040, 1024 x 1024 
pixels^, 30 frames/second), and a high-speed camera for 
measuring velocity distributions (Balsa CADI, 128 x 128 
pixels^, 838 frames/second). For some experiments, we 
obtained high spatial and temporal resolution by using 
the Pulnix-1040 in combination with a stroboscopic LED 
array, with the strobe synchronized so that one light 
pulse occurs at the end of the exposure of one frame 
and the next pulse occurs at the beginning of the subse- 
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quent frame [j, l7| . The most recent data has been taken 
with a Cooke 1200HS camera, which combines high res- 
olution (1280x1024 pixels^) with a high frame rate (625 
frames/second at maximum resolution, higher for lower 
resolution). The acquired digital images are analyzed 
to determine particle locations by calculating intensity 
weighted centers of bright spots identified in the images. 

When the frame rate of camera is faster than the inter- 
particle collision rate, it is possible to measure particle 
velocities by measuring the displacement of the particle 
from one image to the next. The trajectories determined 
from the images are nearly straight lines between colli- 
sions. Although the frame rate is fast compared to the 
mean time between collisions, the fact that it is finite 
introduces unavoidable systematic errors into the exper- 
imentally determined velocity distributions. Some parti- 
cles will undergo a collision between frames, and the re- 
sulting displacement will not represent the true velocity 
of the particle either before or after the collision. Since 
fast particles are more likely to undergo a collision, the 
high energy tails of the distribution will be decreased. 
The probability that a particle of velocity v will not un- 
dergo a collision during a time interval At is proportional 
to exp{—vAt/lo), where lo is the mean free path. The ef- 
fect of this on the velocity distribution function reported 
here is quite small. In addition to taking velocities out of 
the tails of the distribution, collisions will incorrectly add 
the velocities of those particles elsewhere in the distribu- 
tion. This effect can be minimized by filtering out the 
portion of particle trajectories where coUisions occurQ. 

Molecular Dynamics Simulations 

While many of the quantities of interest can be gath- 
ered directly from the experiment, some, like the vertical 
positions and velocities of the particles, are very difficult 
to measure. These same parameters are easy to measure 
using computer simulations which also allow systematic 
variations of material parameters and can provide useful 
confirmation of the interactions underlying the observed 
effects. 

We have written code based on simplified particle in- 
teractions that have been used in molecular dynamics 
simulations of granular media |23 and have successfully 
reproduced, in considerable detail, many of our experi- 
mental results U 0, The inter-particle interactions 
include an elastic repulsive contact force that is propor- 
tional to the particle overlap, an inelastic normal force 
proportional to the relative normal velocities, and a dis- 
sipative tangential friction force proportional to the rel- 
ative tangential velocities of overlapping particles. For 
computational efficiency, the particles in the simulation 
are softer than the steel or brass particles used in the 
experiment, but the parameters give realistic coefficients 
of restitution. This model does not include a tangen- 



tial restoring force that has been found to be necessary 
in some instances [i^, but may be less important in 
the un-sheared systems we have studied so far. In the 
case of simulations of the smooth plate, interactions with 
the top and bottom boundaries are treated like collisions 
with infinitely massive planes, while for the rough plate, 
the bottom boundary is replaced with spheres of infinite 
mass. Either periodic boundaries or elastic or inelastic 
side walls can be used. 

PHASE DIAGRAM 
Free surface layer 

This simple experimental system displays a rich phase 
diagram that depends on the density of particles, the 
plate vibration amplitude A, frequency v, the gap height 
H, and the coefficient of restitution e, of the particles. 
For modest vibration amplitudes, the effect of gravity is 
enough to keep the layer effectively two-dimensional, and 
the lid does not play a significant role. We will discuss 
this regime first and then describe some new results on 
the effect of the confining gap at higher accelerations. 

At very low densities, the interparticle collisions are 
not frequent enough to overcome friction with the bottom 
plate, and the granular layer does not exhibit any signif- 
icant horizontal motion. For a range of surface coverages 
above about p2D = 25%^^, and a range of frequencies 
and accelerations, a fully fluidized, disordered state is 
observed (see j^l for representative phase diagrams). At 
low vibration amplitudes, the fluctuations in the fluidized 
state include long-lived, high density, low granular tem- 
perature (average kinetic energy of horizontal motion per 
particle) clusters @. These fluctuations occur at peak ac- 
celerations at or below Ig, where the dependence of the 
energy flux from the vibrating plate on the dynamics of 
the granular layer ^ probably reinforces the tendency of 
clusters to form due to dissipative interparticle collisions. 

At still lower accelerations, the fluid becomes unsta- 
ble to the formation of a 'collapse', a hexagonal close 
packed layer of particles at rest on the plate, which co- 
exists with a surrounding granular fluid. As with the 
clustering, the mechanism for the instability can be un- 
derstood by considering the energy injection from the 
plate. In particular, at low vibration amplitudes it is 
possible for a particle to remain at rest on the plate, 
but it is also possible for a particle to bounce periodi- 
cally with sufficient energy input to overcome dissipative 
collisions 0, 13 113 > even when F < 1. More generally 
it has been shown that energy injection that depends on 
the local kinetic energy can lead to instabilities in the flu- 
idized stateli^. The collapse-liquid coexistence, and the 
dependence of the pressure in the liquid state on density, 
vibration amplitude, and plate material has been studied 
in detail by Geminard and Laroche 27]. The transition to 
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the fluidized state upon increasing acceleration has been 
studied in Ref. Hi]. 

Above about 80% coverage a fluctuating but hexago- 
nally ordered state is observed for a range of frequencies 
and vibration ampHtudes. First observed by Pieranski 
and co-workers0[ 0|, this phase appears to be sepa- 
rated from the disordered hquid by a continuous phase 
transition that fits the equihbrium 2D melting scenario 
surprisingly well|3d|. In the absence of a confining lid, 
this phase is only observed at relatively high frequencies. 
This may be due to the fact that, for a given accelera- 
tion, at higher frequencies the maximum height of the 
trajectories is lower and therefore the behavior of the 
granular layer is closer to a true two-dimensional system. 
Some evidence for this interpretation comes from the ob- 
servation that the phase diagrams at low accelerations 
measured for two different ball sizes collapses when the 
frequency axis is scaled by the characteristic frequency 
Vc = (g/cr)^/^, where a is the ball diameter [^. The char- 
acteristic time scale, v^^ , is related to the time it takes 
a ball to fall 1 a due to gravity. In the next section we 
show that by confining the layer with a lid while shaking 
with a high vibration amplitude, many of these compli- 
cating gravity-dependent factors can be eliminated, and 
the resulting behavior bears a remarkable resemblance to 
that observed in equilibrium systems. 

The effect of confinement 

We have previously reported the discovery of a dra- 
matic solid-liquid phase transition that occurs in the 
granular layer when the layer is confined with a lid and 
shaken hardQ- Figure |2K shows an image of the two- 
phase coexistence observed in the experiment. This co- 
existence represents a steady state and persists as long 
as the vibration is maintained. The same behavior is ob- 
served in the simulations, and Fig. [2}3, a rendering of 
the simulation data, shows more clearly the structure of 
the ordered phase. It consists of two layers of a square 
packing, rather than the hexagonal packing normally ex- 
pected for hard spheres. 

This surprising behavior can be understood, at least in 
part, by considering the very similar behavior observed 
in equilibrium hard sphere colloidal suspensions confined 
in a similarly narrow gapjs^- These systems show a com- 
plex phase diagram that depends on the volume fraction 
of the colloidal particles and the ratio of the confining gap 
to the particle diameter [ill lil lil lil lilliil . The phase 
diagram can be explained in terms of entropy maximiza- 
tion. The configuration that gives the particles the most 
room to rattle around, within the constraints imposed 
by the gap, will minimize the free energy of the system. 
For a gap of 1.75 a, it turns out that two layer square 
packing is more efficient than hexagonal and hence has a 
lower free energy. As the gap height and volume fraction 
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FIG. 2: A. Experimental observation of coexistence of liq- 
uid and solid phases in a confined granular layer with p2D = 
N/Nmax = 0.9, H = 1.75 cr, 1/ = 80Hz, A=0.085cr. The image 
is an average over 1 second. B. A three-dimensional rendering 
of particle positions near a solid-liquid interface appearing in 
the computer simulation, for similar parameters. Fluid-like 
particles are colored blue, the top layer solid-like particles 
are transparent and the bottom layer solid-like particles are 
colored red. 



are varied, the arrangement of the balls which gives the 
most efficient packing changes. Figure |31 shows several 
of the phase found in our computer simulations of the 
granular layer with different heights and a fixed volume 
fraction. The ordering observed varies from fluid, single 
layer hexagonal (A), buckling (6), two layer squares {2D), 
and two layer hexagonal (2A). There are coexistence re- 
gions between neighboring phases where, for example, 
two hexagonal layers coexist with a fluid (2 A, /). 

We have mapped out the phase behavior of the gran- 
ular system by varying the reduced gap height, h (where 
h — {H — 1) and H is the gap height in units of ball 
diameter, a) and reduced volume density, pH (ph = 
Na^/{SH) where N is the number of particles and S 
is the area of the plate). The phase diagram for the 
granular system, derived from the simulations is shown 
in Fig. 2J3. To determine the phase diagram, simula- 
tions were performed on a grid of ph and h values which 
were initially spaced at regular intervals of 0.05 in both 
variables. After determining the phase at each {ph,h) 
point, further simulations in the transition regions were 
performed with a finer grid spacing of 0.02 in both ph 
and h. Finally, boundary points were interpolated from 
the simulated data and connected with lines to mark the 
boundaries. At values of ph > 1 we are able to observe 
both the 2D and 2A phases without the presence of the 
coexisting fluid. However, for these large values of ph it is 
difficult to accurately map out a phase diagram as there is 
little room for the spheres to rearrange after being placed 
inside the system, so our phase diagram is limited to val- 
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FIG. 3: Four of the phases encountered as a function of re- 
duced gap height, h in the simulations. For clarity, particles 
in the top half of the cell are colored red and particles in the 
bottom half of the cell are colored blue. A. h = 0.3a: Sin- 
gle layer hexagonal. Ji. h — 0.5a: Buckling and single-layer 
hexagonal coexistence. C. h — 0.7a: 2-Layer square and liq- 
uid coexsitence T). h = O.SSct: 2-Layer hexagonal and liquid 
coexistence. All simulations are for ph = 0.9, A = 0.15cr, and 
i/ = 75 Hz. 




ues of p < 1. For comparison. Fig. 2J\ shows the phase 
diagram for the colloidal system studied in Ref. (s^, 
while the granular phase diagram that we found is pre- 
sented in Fig. 2J3. Overall, the phase diagrams are very 
similar, with a few notable differences. The most dra- 
matic difference is the width of the coexistence regions, 
(2n, f) and (2A, f). This difference can be understood in 
terms of the dissipation present in the granular system. 
The crystal structure that forms, whether of square or 
hexagonal symmetry, is of higher density than the sur- 
rounding fluid and hence has a higher rate of dissipation. 
This higher dissipation causes the crystal particles to cool 
and the resulting pressure imbalance compresses the crys- 
tal. Because of the compression the crystal that forms 
in the granular system has a higher density than in the 
colloidal system, and it therefore requires a much larger 
increase in overall density to fill the system with the crys- 
tal phase. We have directly measured a significant dif- 
ference in granular temperature between the coexisting 
phases in both experiment and computer simulations^. 
Newer experiments have shown that the dynamics of the 
spheres inside the crystal is rather complicated. Fig |S1 
shows a typical trajectory for one particle. The numer- 
ous regions where the trajectory is not varying linearly 
with time, imply that there are frequent non-coUisional 
forces on the particle, such as frictional drag forces with 
the lid or enduring contacts with neighboring particles. 

The similarity between the ordering behavior observed 
in the granular and colloidal layers strongly suggests that 



FIG. 4: A. Phase diagram for the hard sphere colloids, re- 
produced from Ref. |3^ . The areas filled in with dashed lines 
represent coexistence regions of neighboring phases. B. Phase 
diagram for the granular system. Here, A = O.lScr, u = 75Hz, 
N = 3000 particles. 
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FIG. 5: Trajectory of single particle in the top layer of the 
two square layer solid phase. The experiment was performed 
using the Cooke 1200HS camera at 4000 frames/second. 
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very similar physics underlies both systems. Unfortu- 
nately, quantitative descriptions of the equilibrium sys- 
tem start with an analysis of the free energy (which, for 
hard spheres, depends only on the entropy). There is 
no recipe for adapting this methodology to strongly non- 
equilibrium granular dynamics, and a useful analysis of 
this phase behavior will have to take into account the role 
of the external driving and energy dissipation. Figure El 
shows the frequency dependence, in both the experiment 
and the simulations, of the vibration amplitude neces- 
sary to nucleate a two square layer crystal. (Note that 
we have plotted the critical value of A, the amplitude 
of the oscillation of the plate position, rather than F.) 
When the vibration amplitude is first increased above 
the critical amplitude, one or more disconnected crys- 
tals quickly nucleate from the homogeneous liquid, far 
from the boundaries of the plate. The crystals grow and 
eventually merge until a steady two-phase coexistence is 
reached. The low frequency increase in the critical 
amplitude occurs because the plate accelerations are not 
large compared to the Earth's gravity, and the particles 
have barely enough energy to reach the second layer. At 
high frequencies, the frequency dependence is very weak 
for steel particles. We have speculated that the formation 
of the solid phase upon increasing vibration amplitude is 
a consequence of layer compression 0. For brass parti- 
cles, with a lower coefficient of restitution e, the critical 
amplitude is greater and there is a slow rise in the criti- 
cal amplitude at higher frequencies. In the simulations, 
we observe the same frequency dependence of the criti- 
cal amplitude and similar dependence on the inelasticity. 
These results show that the nonequilibrium effects due to 
the forcing and the dissipation are an important compo- 
nent of the phase diagram, despite the obvious similarity 
to that observed in the equilibrium colloidal system. 

The phase behavior of highly confined vibrated granu- 
lar particles shows many parallels with the phase behav- 
ior of similarly confined colloidal particles in equilibrium. 
This close correspondence suggests that similar driving 
mechanisms are at work, and that a predictive theory of 
the phase diagram for this system may be derivable from 
the equilibrium framework, with suitable modifications 
to account for the observed dependence on driving and 
dissipation. 



STATISTICAL CHARACTERIZATION OF THE 
GRANULAR GAS 

Introduction 

While developing a kinetic theory of granular flu- 
ids analogous to the kinetic theory of molecular fluids 
at equilibrium 
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has proven to be v ery dif- 
ficult, there has been considerable progress, [iol |4]| 
For example, the hydrodynamic transport coefficents for 
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FIG. 6: The frequency dependence of the critical amplitude 
necessary to nucleate the solid phase. For all cases, p2D = 
0.85, H = 1.75(T. Black squares: steel. Red circles: Brass. 
Green and Blue triangles: MD Simulations 



granular gases have been calculated for some situations 
jU IS 0|. The approach taken to develop the hy- 
drodynamic description of granular gases often begins 
with the solution of the Boltzmann equation in the spa- 
tially homogeneous state, where all gradients are zero. 
In the case of undriven granular fluids, this state is the 
somewhat artificial 'homogeneous cooling state' which 
continually loses energy. The Boltzmann equation for 
this state was first solved in the work by Goldshtein 
and Shapiro ji^- Some recent analytical work and 
Monte Carlo simulations ji^ have focused on the non- 
equilibrium steady state obtained when the energy sup- 
plied by spatially uniform random forcing is balanced by 
collisions. In addition, in a more realistic approach to 
our experimental system, the boundary conditions for a 
vertically vibrated granular fluid have also been studied 
[4^ Isol Isij l and the normal fluid state of the vertically 
vibrated granular layer has been determined, to first or- 
der in the gradients [s^l- Non-Maxwellian velocity dis- 
tributions and velocity correlations have been observed 
in these studies and also in a number of experiments and 
simulations Models incor- 

porating fluctuations in the rate of energy input from 
the driving to the kinetic energy of the particles, either 
through random coUisional transfer[lj, [l^ or through a 
couphng to the fluctuating particle velocities [i^, have 
succeeded in modelling some aspects of the observed dis- 
tributions. Nevertheless, the forcing in the experiments 
is sufficiently different from that of the theoretical models 
and a detailed comparison is difhcult. The effect of corre- 
lations on the transport coefficients is still being studied 
152. I53I I54L I55I l56l ISTl. Below we review our results 
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on the statistical characterization of a vertically driven 
granular monolayer and present measurements of velocity 
autocorrelation functions. 



Pair Correlation Functions 

Correlations in particle positions are most directly 
measured by the pair correlation function, G(r): 



G(r) = 



< p{r')p{r' + r) > 
< p(r') >2 ' 



(1) 



where p is the particle density and the average is taken 
over the spatial variable r'. The correlations of a two- 
dimensional gas of elastic hard disks in equilibrium are 
due only to geometric factors of excluded volume and 
are independent of temperature. This provides a useful 
reference point for the measurements in the experimental 
system. 

We have observed three distinct types of behavior 
of the pair correlation function. At low vibration am- 
plitudes where strong clustering is observed, G{r) has 
deeper modulations and a higher value at contact than a 
system of equilibrium hard disks at the same density 0. 
The relatively long-lived high density, low granular tem- 
perature fluctuations described in Sec. b are presumably 
responsible for the strong spatial correlations. 

At plate accelerations considerably above Ig in the sys- 
tem without a hd, G(r) loses all of its structure as the 
layer goes fully three-dimensional. The change to 3D be- 
havior can be seen by the increase in the value of G(r) 
for r < 1 . (The measured correlation function includes 
only horizontal particle separations, so that colliding par- 
ticles at different heights appear to interpenetrate.) This 
change can affect the dynamics in several ways. The 
effective density is decreased so that excluded volume ef- 
fects are less important, the inter-particle collisions can 
occur at angles closer to vertical, affecting transfer of en- 
ergy and momentum from the vertical direction to the 
horizontal, and the change in dimensionality itself can 
have important consequences. 

In order to separate these effects from the direct con- 
sequences of increasing the kinetic energy of the gas, the 
Plexiglas lid was added to the system, with a gap small 
enough so that the particles cannot pass over top of one 
another, although enough room remains for collisions be- 
tween particles at sufficiently different heights to transfer 
momentum from the vertical to the horizontal direction. 

G{r) measured with the lid on shows that the particle- 
particle correlations persist and become independent 
of r when the system is constrained in the vertical 
direction^. The structure observed in the correlation 
function is indistinguishable from that of an equilibrium 
elastic hard disk gas at the same two-dimensional den- 
sity, indicating that the correlations that exist are due 



to excluded volume effects. From F = 1.50 to F = 3.0, 
the granular temperature changes by more than a fac- 
tor of 2f3|, yet there is no detectable change in the pair 
correlation function. 



Velocity Distributions 

The distribution of horizontal velocities of the gran- 
ular particles is typically non-Maxwellian. As with the 
pair correlation function, we have observed a number of 
different regimes. At low vibration amplitudes where the 
clustering is very strong, the velocity distribution has 
a nearly exponential shape, with a cusp at zero veloc- 
ity and high velocity tails that are nearly straight on 
a log-linear plotQ. The broad distribution presumably 
arises from the fact that the local granular tempera- 
ture, like the local density, is inhomogeneous in the clus- 
tering regime Q, and it has been shown very recently 
that the cusp at zero velocity is likely due to friction 
with the bottom plate [s^. As the acceleration is in- 
creased, the velocity distribution crosses over smoothly 
to a Maxwellian distribution |^. Identical behavior was 
observed in molecular dynamics simulations designed to 
mimic the experimental situation[23|. In addition, those 
simulations showed that, at low accelerations, the aver- 
age kinetic energy (granular temperature) in the vertical 
motion was several times larger than that in the hori- 
zontal. The approach to Maxwellian horizontal velocity 
distributions coincides with a dramatic reduction in the 
anisotropy. This is consistent with the observation that 
non-Maxwellian features in the velocity distribution can 
arise from energy transfer through collisions in strongly 
anisotropic situations|lfil|. The expansion of the layer in 
the vertical direction as the acceleration is increased pre- 
sumably also results in more interparticle collisions rela- 
tive to particle-plate collisions, and this ratio may play 
an important role in determining the form of the velocity 
distributions 01 . 

When the expansion of the layer is eliminated by 
adding the lid, the crossover to Maxwellian distributions 
is not observed^. As with G{r), the velocity distribu- 
tion becomes essentially independent of the acceleration 
for values above roughly 1.5 g. However, whereas G(r) 
becomes essentially indistinguishable from that observed 
in equilibrium systems, the velocity distribution has a 
cusp at zero velocity and high-velocity tails that are sig- 
nificantly broader than Gaussian tails. 

When the smooth bottom plate is replaced by the 
rough plate, the cusp at zero velocity disappears, but the 
broad high-velocity tails remain, at least for moderate 
densities ?]. At very low densities, however, the distribu- 
tions approach a Maxwellian distribution. Measurements 
at low densities (typically below p2D = 0.25) are not pos- 
sible on the flat plate, because the interparticle collisions 
are too infrequent and all horizontal motion ceases. By 
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contrast, every collision with the rough plate injects en- 
ergy directly into the horizontal, as well as vertical, ball 
motion. As a result even a single ball on the plate will ex- 
hibit apparently random changes in direction, much like 
the random forcing of some models of heated granular 
fluids [11 m H ll9j. As described below, the behavior 
of the velocity autocorrelation function suggests that the 
scattering of the ball with the rough plate acts as a vis- 
cous drag. Thus the Maxwellian velocity distributions 
at low densities may result from Langevin-like dynamics, 
with the plate providing both the random forcing and 
the drag. As the density is increased, both the drag and 
the interparticle collisions will affect the velocity distri- 
bution, a situation modeled by Puglisi, et al.[l^. How- 
ever the granular temperature is always anisotropic in the 
quasi-2D layers, so the energy transfer from the vertical 
to the horizontal components of the velocity will affect 
the velocity distribution on the rough plate as well llSl. A 
somewhat similar geometry was investigated in Ref. |59( , 
where a full layer of heavy two-particle chains played the 
role of the rough plate. This provides a dynamic bottom 
layer, in contrast to our rough plate, where the bottom 
layer of balls is glued to the oscillating flat surface. The 
velocity distributions appear to be closer to a Maxwellian 
than those we have reported (compare Fig. 5 of Ref 
with Fig. 1 of Ref. |E^), but there are a number of 
differences between the experiments, including the rela- 
tive sizes of the balls in the two layers and the material 
properties of the particles, so that it is not clear if the 
quantitative differences are due to the motion of the bot- 
tom layer or to other effects. 



Velocity autocorrelation function 

As described above, continuum models of granular 
flows need to account for the presence of correlations, 
which can be seen directly in the self-diffusion of individ- 
ual particles. We have measured the normalized velocity 
autocorrelation function. 



A{t) 



> 



(2) 



where Vi is one component of the horizontal velocity and 
the brackets represent the time average. While in prin- 
ciple the hexagonal symmetry of the rough plate could 
produce an asymmetry in the autocorrelation function, 
none is observed, so the data shown is the average of 
the result for each component separately. Fig. [7^ shows 
A{t) measured on the rough plate at different densities 
at the same plate acceleration. As expected, the auto- 
correlation function decays more slowly as the density 
decreases because of the longer time between interparti- 
cle collisions. At the lowest densities, however, the decay 
rate becomes almost independent of the density, suggest- 



ing that scattering from the rough plate is the dominant 
mechanism of momentum decay at those densities. 

The autocorrelation functions for brass and steel 
spheres at a density of p2D ~ 0.6 are compared in Fig. [3d. 
We used a slightly higher plate acceleration for steel than 
for the more dissipative brass (F = 1.57 vs. F = 1.39) to 
produce the same horizontal granular temperature. The 
autocorrelation functions show an approximately expo- 
nential decay at short times, as it can be seen more 
clearly in Fig. Ct, when plotted on a log- linear plot. 
The decay rate is about a factor of two faster than that 
calculated from kinetic theory for a liquid at the same 
density and temperature, presumably due to the added 
effect of the scattering with the rough plate discussed 
above. At longer times there is clearly upward curvature 
on the log-linear plot. This slow temporal decay is likely 
related to the equal time spatial correlations reported 
previouslyQ and described below. (Although long time 
tails in the velocity autocorrelation function are present 
in dense elastic fluids as well'H^, those effects are smaller 
than what we observe in the granular fluid ) . The slower 
decay is more dramatic for the brass, probably because 
the spatial velocity correlations, which are a consequence 
of the forcing and dissipation, are stronger, but we have 
not confirmed this. 

The diffusion coefficient is directly related to the inte- 
gral of the velocity autocorrelation function, but it can 
be measured most readily from a plot such as Fig. [3i, 
which shows the mean squared particle displacement di- 
vided by time for the brass and steel spheres. For normal 
diffusive behavior, this curve will reach a constant value 
(proportional to the diffusion coefficient) at long times. 
Both brass and steel clearly exhibit normal diffusion, but 
the diffusion coefficient for brass is higher than that of 
steel. In the absence of correlations, the diffusion coeffi- 
cient given by hard-disk kinetic theory depends only on 
the mean free path, which depends only on the density, 
and the granular temperature. Thus the enhanced dif- 
fusion in brass is a direct demonstration of the effect of 
forcing and dissipation on transport properties in granu- 
lar gasses. 



Spatial velocity correlations 

The assumption of 'molecular chaos', that the veloci- 
ties of colliding particles are uncorrelated, is a crucial ap- 
proximation normally used to solve the Boltzmann equa- 
tion and to calculate other fundamental quantities in the 
kinetic theory of liquids. The presence of spatial velocity 
correlations can significantly affect coUisional transport 
of mass, energy, and momentum, [islfiil 

The longitudinal and transverse velocity correlations. 



9 




i:(ms) -[(ms) 




FIG. 7: Autocorrelation and self-diffusion on the rough plate A. A{t) at different values of p2D at F = 1.5 B. A{t) 
for steel and brass at p2D ~ 0.6 and F = 1.39 (steel) and F = 1.57 (brass). The vibration amplitude was chosen so that the 
granular temperature was the same in both cases. C. Same data as in B on a log-linear plot. D. Mean square displacement 
divided by time, for the same conditions as B. All data were collected for iy—60 Hz. 
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C|| and Cj_, respectively, are calculated by 

N,. 

where the sum runs over the iV^ pairs of particles sepa- 
rated by a distance r, wj' is the projection of along the 
line connecting the centers of particles i and j, and is 
the projection perpendicular to that line. (Only the hor- 
izontal components of the velocity are measured). Mea- 
surements on the smooth plate revealed surprising strong 
negative values for C[|(r), indicating that the particles in 
the granular gas exhibit a tendency to move apart from 
one another|3- On the smooth plate, there is always 
more energy in the vertical components of the velocity 
than in the horizontal components, because energy is in- 
jected exclusively into the vertical motion. This energy 
gets transferred to the horizontal components during col- 
lisions, so that the particles have a tendency to 'fly apart' 
from collisions, resulting in the observed anti-correlated 
velocities. 

Measurements on the rough plate, by contrast, show 
C|| _L(r) everywhere positive, with an exponential decay 
that is approximately independent of density. A model 
system with random forcing also shows positive correla- 
tions, but with a slower decay that depends on the den- 
sity of the granular gas 1^ 0, 0, |^. This difference 
may be due to the scattering of the balls by the rough 
plate. These collisions do not conserve momentum, and 
are therefore probably more effective than ball-ball colli- 
sions at dissipating the velocity correlations which arise 
from the effects of the random forcing. The dramatic 
change in the velocity correlations when the bottom sur- 
face is changed provides a direct demonstration that the 
kinetics of the granular gas are very sensitive to the forc- 
ing mechanism. 

CONCLUSION 

The surprisingly rich behavior of the simple system 
described in this paper is in many ways emblematic of 
the opportunities and challenges presented by granular 
flows more generally. There are situations where the sta- 
tistical mechanics of thin vibrated layers is essentially 
similar to equilibrium systems, such as the existence of 
nearly Maxwellian velocity distributions at high vibra- 
tion amplitudes in the absence of confinement, or at low 
densities on the rough plate. On the other hand, there 
are situations where the driving and the dissipation pro- 
duce strongly non-equilibrium behavior, such as strongly 
non-Maxwellian velocity distributions and significant ve- 
locity correlations, evidence for violations of molecular 
chaos. The observation of a complex phase diagram that 
is clearly similar to a well understood equilibrium sys- 
tem suggests that entropy maximization, the mechanism 



that drives the transitions in equilibrium, can operate 
far from equilibrium. Adapting that formalism to ac- 
count for the observed effects of forcing and dissipation 
remains an open challenge. 
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